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We compare the three most commonly used scanning probe techniques to obtain a reliable value of the work 
function in graphene domains of different thickness. The surface potential (SP) of graphene is directly 
measured in Hall bar geometry via a combination of electrical functional microscopy and spectroscopy 
techniques, which enables calibrated work function measurements of graphene domains in ambient 
conditions with values 0 1LG —4.55 ± 0.02 eV and <P 2 lg ~ 4.44 ± 0.02 eV for single- and bi-layer, 
respectively. We demonstrate that frequency-modulated Kelvin probe force microscopy (FM-KPFM) 
provides more accurate measurement of the SP than amplitude-modulated (AM) -KPFM. The discrepancy 
between experimental results obtained by different techniques is discussed. In addition, we use FM-KPFM 
for contactless measurements of the specific components of the device resistance. We show a strong 
non-Ohmic behavior of the electrode-graphene contact resistance and extract the graphene channel 
resistivity. 



Mapping of the local electronic properties of graphene is necessary for control of growth parameters and 
for understanding device functionality. The accurate quantification of the measured values is essential in 
order for the properties of graphene to be reliably understood and compared. 
The growth of graphene by the sublimation of Si from SiC is arguably the most advanced method for producing 
continuous, homogeneous large area graphene 1 . Control over the layer thickness has been demonstrated 2 , with 
sublimation being the method of choice for device manufacture, where a continuous, large area of single-layer 
(1LG) of graphene is required, i.e. nanoelectronics, sensing, THz applications, etc. Due to advances in the sample 
growth, it is generally possible to achieve homogenous single layer coverage over large areas (i.e. —95% 1LG 
coverage for the sample presented in this work). However, even small inclusions of bi-layer graphene (2LG) leads 
to redistribution of carriers, inhomogeneous screening effects, and the corresponding nanoscale changes in the 
surface potential (SP) and work function (&). Unambiguous determination of the layer thickness of epitaxially 
grown graphene using atomic force microscopy ( AFM) is particularly challenging due to the stepped nature of the 
SiC substrate coupled with growth of graphene, which often nucleates at step edges 1 . 

Scanning measurement techniques, such as Kelvin probe force microscopy (KPFM), are widely used for 
mapping the SP of graphene as well as identification of graphene layers. For example, KPFM has recently been 
used to distinguish between areas of 1LG, 2LG, few layer graphene (FLG) and the buffer or interfacial layer (OLG) 
and for validation of optical quality control methods for graphene 3,4 . However, KPFM does not generally provide 
reliably comparable values for differences in SP between layers with a wide variety of AV C pd 1 ' 2LG values previously 
reported for 1-2LG. For example, for epitaxial graphene on SiC, Filleter et aV reported a Ay CPD 12LG =135 mV in 
vacuum, whereas Burnett et at. 5 obtained a Ay CPD 12LG = 25 mV in air. On the other hand, in case of exfoliated 
graphene on Si0 2 , Yu et at. 6 reported a Ay CPD 12LG = 120 mV after accounting for environmental effects by 
measuring in ambient atmosphere and dry nitrogen, whereas Ziegler et al. 7 reported a smaller value of Ay CPD 12LG 
= 68 mV in ambient conditions. 

The effects of substrate on the charge transfer to graphene and subsequent change of the SP have been 
discussed in depth 7 " 10 . Additionally, change in the charge carrier concentration, whether it is intentional, 
by electrostatic or photochemical gates 1112 , or incidental, as by uncontrolled adsorbates, modifies the 
measured Ay CPD 12LG values (see ref. 6). For example, specific atmospheric gating can modify the 
AV C pd 1 ' 2LG of epitaxial graphene from 0 to 100 mV on changing the environment from vacuum or pure 
nitrogen to >1 ppm N0 2 in nitrogen mixture 13 . Moreover, atmospheric humidity gating has been shown 
to increase the Ay CPD 12LG values 1415 . While the reported discrepancy in the published values of Ay CPD 12LG 
can be partly attributed to different substrate and environmental gating, here we primarily address the 
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measurement methodology and the accuracy of KPFM measure- 
ment technique applied to graphene domains. 

We analyze the ability to obtain quantified, comparable and 
accurate results of single-pass frequency-modulated (FM)-KPFM, 
conventional dual-pass amplitude-modulated (AM) -KPFM and 
electrostatic force spectroscopy (EFS) by performing measurements 
on a graphene Hall bar device after SP calibration of the AFM probe 
against gold electrodes. In contrast to many experimental studies, we 
aim to investigate graphene devices in standard ambient conditions 
(rather than in vacuum or specific gas atmosphere), as such condi- 
tions are the most representative both for general research and indus- 
trial lines. We find that conventional AM-KPFM, being a force 
sensitive technique, suffers from a spatial averaging effect of the SP 
due to a significant contribution of the cantilever base and cone to the 
capacitive coupling, which reduces Ay CPD 12LG and leads to incorrect 
values of SP measured on a biased device. In contrast, FM-KPFM is 
sensitive to the force gradient and measures the SP of the area directly 



under the probe apex, demonstrating improved spatial resolution 
and absence of averaging effects. We conclude that FM techniques, 
such as FM-KPFM and EFS, provide more accurate measurement of 
the SP than AM-KPFM. Using calibrated FM-KPFM, we perform 
precise work function measurements of 1LG and 2LG, being @i LG = 
4.55 ± 0.02 eV and 0 2L g = 4.44 ± 0.02 eV, respectively, for the 
sample studied here. We also perform contactless measurements of 
the resistance of the graphene channel and two separate electrode - 
graphene lead contacts. 

We demonstrate that the experimental approach presented here 
can be successfully used for standardization of measurements of 
the work function and obtaining reliable quantitative parameters 
not only in graphene but in many other electronic materials, i.e. 
semiconductors, photovoltaic, etc. Representing a surface state of 
a material rather than its bulk property, the work function, in 
graphene in particular, can be strongly affected by environmental 
conditions. To assure accuracy of measurements, recalibration 
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Figure 1 | Schematic diagrams of the experimental techniques, (a) AM-KPFM; (b) FM-KPFM; topography of the graphene Hall bar is superimposed 
with SP maps on a 3D image. Plots show characteristic profiles, i.e. SP on top and topography on bottom along the horizontal line in the center 
of the image (not shown), (c) Typical parabolic change of the cantilever phase shift measured by EFS during DC voltage sweep at a fixed point on 1LG. 
(d) UPS data showing the work function of gold, @ Au = 4.82 eV, for four different samples. 



SCIENTIFIC REPORTS | 3 : 2597 | DOI: 1 0.1 038/srep02597 



2 



should be done upon any significant change of ambient (e.g. 
humidity). 

Results 

Surface potential measurement techniques. Surface potential maps 
of a sample can be obtained using KPFM, which measures the 
strength of the electrostatic forces between a conductive probe and 
the sample 5 . There are different methods of detecting electrostatic 
forces, namely: AM-KPFM, which responds to the electrostatic force 
at a set frequency of probe oscillation (Figure la); and FM-KPFM, 
which responds to the electrostatic force, while maintaining constant 
amplitude of cantilever oscillation (Figure lb). As we show below, the 
choice of the measurement technique significantly affects the 
accuracy of surface potential measurements on micrometer scale 
graphene. The schematic diagrams of the used techniques are 
shown in Figure 1; all techniques are discussed in detail in the 
Method section. 

AM-KPFM: experimental results. Figure 2a shows a topography 
map of the graphene device. SiC step terraces are clearly visible 
running at a —60° angle to the channel. Gold contacts are seen at 
the left and right hand sides of the image. The image reveals that it is 
generally rather difficult to determine the graphene layer thickness 
from the topography maps. The surface potential of the electrically 
grounded device was mapped using AM-KPFM in ambient 
environment (Figure 2b). Bright areas of 2LG are clearly visible on 
the 1LG background, whereas darker regions correspond to etched 
SiC. The value of the AV C pd 1 ' 2LG = 50 mV is consistently measured 
over all areas of the sample (Figure 2c), which is comparable to 
previously published results on similar samples 8 . The SP dip which 
can be seen to the right of the 2LG step in Figure 2c (as well as in 
Figure 3b and Figure 4b) is attributed to a small patch of resist 



residue, which is clearly observed in the topography map (circled 
in red in Figure 2a). 

Further to this, we study the surface potential of a biased graphene 
device. Bias voltages of V c h = 0, ±0.5, ±1, ±1.5 and ±2 V were 
applied to the left gold electrode and SP maps of the device were 
obtained in AM-KPFM mode. Figure 2d shows the plotted SP values 
along the marked line (Figure 2b) going through the center of the 
channel and connecting the gold leads. The raw data is plotted in 
Figures 2c and 2d, i.e. no calibration of the probe's work function has 
been performed here. As a result, the measured SP values in Figure 2d 
are not centered at 0 V. A significant discrepancy between applied 
and measured voltages is observed using AM-KPFM, i.e. the total 
difference in surface potential values measured on the left gold elec- 
trode, when biased with V ch = +2 and —2 V, is only —2.9 V, i.e. 
27.6% less than the expected 4 V. After taking into account the work 
function of the probe, the values of AV C pd between the biased gold 
contacts are still smaller than expected. This discrepancy in applied 
and measured voltages can be explained by the spatial averaging of 
AM-KPFM due to the long-range nature of the electrostatic forces 
acting on the probe and leading to substantial contributions from the 
probe cone and the base 16 . These parasitic contributions can affect 
the measured SP, as the area under the cantilever may not be directly 
over the gold leads, but instead averaging the SP over the graphene 
device leading to a lower total value. For a given device geometry and 
using AM-KPFM, it is expected that scanning across the channel 
might somewhat decrease parasitic capacitive coupling between the 
cantilever and the gold electrodes 16 . However, bearing in mind that 
the size of the cantilever (200 X 30 um 2 ) is considerably larger than 
the device channel (50 X 5 um 2 ) and intricate device electrode geo- 
metry (up to 6 electrodes and bonding pads of a complex shape), 
some coupling between the cantilever and the gold electrodes and 
bonding pads is unavoidable in any scanning direction. Moreover, 
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Figure 2 | Topography and surface potential mapping with AM-KPFM. (a) Topography map of the device showing a double-cross Hall bar, gold 
electrodes and etched areas (bare SiC) which define the channel. Red circle denotes a small fraction of the resist residue, which is seen as a dip on all SP 
profiles, (b) AM-KPFM surface potential map of the grounded Hall bar device, (c) Plot of the surface potential between areas of 1LG and 2LG within the 
channel along the dashed line shown in the inset. Inset shows the magnified area of the AM-KPFM surface potential map framed in (b). (d) Plot of the 
surface potential for the biased device measured between gold leads through the center of the channel along the line depicted in (b), the left gold lead is 
biased at between +2 and —2 V and the right gold lead is grounded. 
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any inhomogeneity of the device itself (i.e. the presence of graphene 
domains of different thickness) will contribute to the averaging effect 
of AM-KPFM technique. It should be noted that the maintaining of 
absolute uniformity of graphene thickness over the length of 200 um 
remains challenging. 

FM-KPFM: experimental results. Surface potential mapping has 
been further carried out using FM-KPFM on the same device. 
Figure 3a shows the potential map of the grounded device. Areas 
of 1LG and 2LG are sharply outlined and better defined compared to 
the measurements taken with AM-KPFM. Values of AV C pd' 2LG = 
150 mV are recorded as shown in Figure 3b. This value is consistent 
over the device and significantly larger than Ay CPD 12LG obtained 
with AM-KPFM. The larger Ay CPD 12LG values can be accounted 
for by considering the measurement technique, which uses the 
force gradient rather than the force and also leads to improved 
spatial resolution of FM-KPFM compared to AM-KPFM. Figure 3c 
shows a line profile of the surface potential measured along the center 
of the channel with V ch = 0, ±0.5, ± 1, ± 1.5 and ±2 V. The change 
in surface potential values measured on the left gold lead when biased 
with V ch = +2 and -2 V is now -4.18 V, i.e. 4.5% larger than the 
expected 4 V, suggesting that this technique provides improved SP 
measurements even over relatively small structures with a size of 
several micrometers. This result is in a very good agreement with 
recent finding 17 , where a negligible averaging effect was demon- 
strated for graphene samples using FM-KPFM mode. 

Electrostatic force microscopy and spectroscopy: experimental 
results. Figure 4a shows an EFM phase map of the grounded 
device. Due to the high spatial resolution, the edges of 2LG 
domains are sharp and well defined. Figure 4b shows the recorded 
SP values obtained from EFS measurement points over the area of 
1LG and 2LG. As this is a spectroscopy rather than a mapping 
technique, values may be slightly affected by the exact position. In 
this instance, Ay CPD 12LG = 110 mV is in a reasonable agreement 
with the results obtained by FM-KPFM, this being expected as both 
techniques are sensitive to the force gradient. The discrepancy can be 
attributed to only a few EFS experimental points obtained on the 
small isolated 2LG domain at the center of the channel, whereas a 
significantly larger number of points were measured with FM- 
KPFM. Further improvement of EFS method and better agreement 



with FM-KPFM can be achieved by decreasing the step size between 
measurement points. Results of measurements of 200 EFS 
spectroscopy points taken along the center of the channel between 
the two gold contacts with the left contact biased at V ch = 0, ± 1 and 
±2 V are shown in Figure 4c. 

Even the most accurate FM-KPFM and EFS techniques provide a 
non-zero reading of the surface potential on the grounded electrode, 
i.e. V CPD = -365 mV for FM-KPFM (Figure 3c) and V CPD = 

— 723 mV for EFS (Figure 4c). This discrepancy is the result of a 
work function difference between the gold and PFQNE-AL probe. 
Further to this, we account for the resulting work function difference 
by subtracting the V ch obtained from the grounded right contact 
from the experimental value of the V C pd> i>e> AV = V CPD (V ch ) — 
Vcpd(0). The procedure was performed using results of all three 
experimental techniques for the range of applied V chi providing 
Ay for the left gold electrode (Figure 5). The measured potential 
drop is typically 27.6% lower than the actual V ch for AM-KPFM, 
whereas it is 4.4% and 7.8% higher for FM-KPFM and EFS, respect- 
ively. The lower Admeasurements are consistent with spatial aver- 
aging, as the relatively large base of the cantilever weakly interacts 
with the device channel and the right contact 16 , both of which are at a 
lower V C pd than the left contact, as was discussed above. The higher 
Admeasurement with FM-KPFM could be a result of an overestima- 
tion of the SP due to a relatively large excitation voltage of V AC = 8 V, 
whereas the discrepancy with EFS rises from un-optimized fitting 
parameters. 

Work function calibration. We employed the use of force gradient 
techniques to provide accurate measurements of work function of 
1LG and 2LG. Initially, work function of the PFQNE-AL probe was 
calibrated against the work function of the gold leads: ^ pro be ~ ^a u 
+ eAVcp D , where Vcpd was measured on the grounded gold 
electrodes. Ultraviolet photoelectron spectroscopy (UPS) measure- 
ments were carried out on four separate samples of gold deposited by 
e-beam evaporation under the same conditions as the deposition of 
the gold electrodes. The spectra were acquired with voltage of 

— 19.04 V applied to the sample. The Fermi edge was centered at 
0 eV by measuring the offset from a high resolution Fermi edge 
spectrum of the silver calibration sample. The offset was used to 
correct the energy scale for all four Au spectra (Figure Id). Using 
the indicated energies obtained from the spectra from each area of 
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Figure 3 | Surface potential mapping with FM-KPFM. (a) FM-KPFM surface potential map of the grounded Hall bar device, (b) Plot of the surface 
potential between areas of 1LG and 2LG within the channel along the dashed line shown in the inset. Inset shows the magnified area of the FM-KPFM 
surface potential map framed in (a), (c) Plot of the surface potential measured between gold leads through the center of the channel along the line depicted 
in (a), the left gold lead is biased at V c h between +2 and —2 V and the right gold lead is grounded. 
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Figure 4 | EFM mapping and surface potential line profiles with EFS. (a) EFM phase map of the grounded Hall bar device, (b) EFS plot of the 
surface potential between areas of 1LG and 2LG within the channel along the dashed line shown in the inset. Inset shows the magnified area of the EFM 
phase map framed in (a), (c) Plot of the surface potential measured by EFS between gold leads through the center of the channel along the line depicted in 
(a), the left gold lead is biased at V c h between +2 and —2 V and the right gold lead is grounded. 



the samples, the work function was calculated. The difference in 
energy between the Fermi edge measured on a silver calibration 
sample and the cut off (x) is given by x = h — 0, where the energy 
of the incident photon is h = 21.22 eV. The cut off was obtained by 
fitting a line to the relevant part of each spectrum, determining its 
gradient and the point at which it crosses the energy-axis. The 
absolute value of x = 16.40 eV was established, thus measuring the 
work function of all four gold samples as 0 Au = 4.82 eV. UPS 
characterization performed in ultra-high vacuum (UHV) would 
include all irreversibly bound adsorbates (chemisorbed oxygen and 
physisorbed hydrocarbon) attached to the surface prior to the 
measurements. Some uncertainty may arise from the reversible 
adsorption of species, such as water, which may form surface 
dipoles leading to a small change of the 0 Au on transferring the 
sample from UHV to ambient environmental conditions. Owing to 
its relatively high standard electrode potential of 1.52 eV 1819 , gold is 
very stable in air, being less prone to both oxidation and formation of 
a submonolayer of water in ambient conditions (both these factors 
could potentially affect 0 Au ) than most other good conductors. 
Following Ref. 20, we can estimate that 0 Au may decrease by —3% 
(top estimation) as the relative humidity changes from 0 to 40%. Use 
of a more stable metal would be ideal for tip calibration, however gold 
remains one of the best electrode materials. 
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Figure 5 | Direct comparison of the surface potential measurements. 

Normalized SP values as measured by AM-KPFM, FM-KPFM and EFS 
techniques on the left gold electrode in dependence on the voltage applied 
to the same gold electrode. 



Using the measured work function of gold, we calculated the work 
function of the probe to be 0 pro be = 4.09 eV using EFS (Figure lc). 
Then, the work function of 1LG and 2LG were determined: 0 sam pie ~ 
0 probe ~ e^VcPD using the measured values of the SP extracted from 
the line profiles of the potential maps, i.e. Vcpd = — 454 and 
— 344 mV for 1LG and 2LG, respectively, see Figure 4b. This defines 
work functions of 0 1LG ~ 4.55 ± 0.02 eV and 0 2L g ~ 4.44 ± 
0.02 eV. These work function values are within the range of prev- 
iously reported results of 4.41-4.57 eV for 1LG measured with FM- 
KPFM 6,21 . It should be noted that the work function of graphene 
depends on the carrier density and is, therefore exceptionally sens- 
itive to substrate and environmental gating due to its two-dimen- 
sional nature. For example, the published values were reported to 
change with varying lab ambient, i.e. the change of work function and 
SP due to adsorbates being <50 meV 6 and —130 mV 21 , respectively. 

Contactless resistance measurements. High accuracy of FM-KPFM 
technique provides an excellent contactless method for measuring 
the electrode-graphene contact resistance with no need for 
specifically patterned electrodes 6 , which is typically used with the 
transmission line method. Using experimental results shown in 
Figure 3c (i.e. line profiles of V CPD at V ch = ±2 V), contact and 
channel resistance can easily be deduced by normalizing these line 
profiles [V C pD(Vch) ~ VcPD(0)]IV c h = & V/V ch as shown in Figure 6a. 
This procedure accounts for any intrinsic V CPD changes, i.e. 
variations in the work function of features, such as 1LG, 2LG and 
gold. The resulting normalized line profile is solely a consequence of 
the potential drop at electrode-graphene contacts and along the 
graphene channel due to changes in the resistance. Dependence of 
the normalized voltage drop AV/V ch as measured across the left 
(right) contacts and graphene as well as across the graphene 
channel (i.e. points 1-2, 3-4 and 2-3, respectively, in Figure 6a) 
are plotted in Figure 6b. While the voltage drop within the 
graphene channel is constant for all applied V ch , this value changes 
linearly on electrode-graphene contacts. Careful inspec- 
tion of the electrode-graphene potential drop for both contacts 
reveals a clear V ch dependence. Focusing on the left contact (points 
1-2), relative change of the voltage on electrode-graphene channel is 
Ay = 0.55 and -0.91 V for V ch = +2 V and -2 V, respectively. 
However, at the right contact, the Ay = 0.83 and -0.52 V for V ch = 
+2 V and —2 V, respectively. From potential drop and I-V ch 
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Figure 6 | Contactless resistance measurements with FM-KPFM. 

(a) Normalized surface potential line profiles. Experimental values are 
obtained by FM-KPFM along the dashed line in Figure 3a. (b) A VI V ch and 
(c) resistance measurements of left contact, right contact and across the 
graphene channel, i.e. points 1-2, points 3-4 and points 2-3, respectively, 
in (a). Dashed lines are guides for the eye. Inset in (b) shows the 
dependence of the total current (I) through the circuit on the bias voltage 
(V ch ). 



(transport) measurements (Figure 6b main panel and inset, 
respectively), the contactless resistance can be determined as 
AV/I(V ch ). Figure 6c shows the contactless resistance 
measurements of the graphene channel (R ch ) only and electrode- 
graphene contacts (R CO nt) separately for the left and right contacts. 
The resistance of the graphene channel R ch — 33 kQ is constant over 
the range of applied voltages, i.e. independent on the V ch . The 
corresponding resistivity value is p ch ~ 2.7 X 10" 6 Ohm cm. On 
the other hand, the contacts exhibit a significant change of AR ie f t 
com ~ — 17.5 kQ and AR right cont ~ 17.0 kQ as V ch changes from 
—2 to +2 V for left and right contacts, respectively. These 
measurements show that for this particular device, R cont is 
dependent on the V ch revealing a non-Ohmic behavior. It should 
be noted that the edges of both gold electrodes overlap onto 2LG 
islands (Figure 3b), which exhibit a lower work function. The lower 
@2lg is a result of a high carrier density of the 2LG system, which 



could affect the flow of charge and, thus the contact resistance. 
However, it is rather difficult to give a quantitative measure of this 
effect from our experiments. 

Discussion 

Using epitaxial graphene Hall bars with gold electrodes, we have 
demonstrated significant differences in accuracy and resolution 
between AM-KPFM and FM-KPFM techniques in determining the 
difference in surface potential between 1LG and 2LG. Values of 
Ay CPD 12LG measured with FM-KPFM demonstrate a threefold 
increase as compared to AM-KPFM. While AM-KPFM measures 
the electrostatic force on the cantilever, FM-KPFM is sensitive to 
the force gradient. Thus, AM-KPFM gives a weighted average of 
the signal, including contributions from the surface under the probe 
cone and cantilever. Sensitivity to the shorter range force gradient 
characteristic for FM-KPFM leads to spatially- confined contribu- 
tions, which arise only from the probe apex, thereby reducing the 
spatial averaging of measured work functions observed in AM- 
KPFM. We experimentally demonstrate that FM-KPFM and conse- 
quently EFS have a greater degree of spatial resolution (<20 nm) 
than AM-KPFM. Improvement in spatial resolution is clear on com- 
paring the sharpness of the potential maps obtained. Moreover, we 
show that use of FM-KPFM and a calibrated probe provide a simple 
and straightforward method of obtaining an accurate measure of the 
work function of 1LG and 2LG. Accuracy of measurements is pro- 
vided by initial calibration of the KPFM signal against the gold elec- 
trodes, whose work function was measured independently by UPS, 
however, keeping in mind that the work function of gold may some- 
what change on transferring the sample to ambient due to adsorption 
of reversible species as discussed above. This improvement in mea- 
surement technique enables greater accuracy in determination of the 
work function of 1LG and 2LG, with values of @ 1L g ~ 4.55 ± 
0.02 eV and @ 2 lg ~ 4.44 ± 0.02 eV, respectively, as valid for the 
particular sample studied here and in specified ambient conditions. 
Thus, the experimental procedure outlined in this paper can be 
implemented as a generic route for standardization of work function 
measurements applicable not only to graphene but to a wide class of 
nanoscale materials and devices, for example in light emission and 
solar devices, displays, etc. 

FM-KPFM was also used to investigate: i) the contact resistance 
between the gold electrode and graphene, revealing a non-Ohmic 
behavior, and ii) the resistance of the graphene channel showing 
Ohmic behavior with R ch ~ 33 kQ and p ch ~ 2.7 X 10" 6 Ohm 
cm. This simple contactless method can be used to investigate the 
specific components of the total resistance, without fabricating 
devices for the transmission line method. 

Our results unambiguously demonstrate that the measured values 
of the SP of single- and bi-layer graphene largely depend on the 
accuracy of the technique. However, we would like to stress that 
the obtained absolute values are not by any means the fundamental 
parameters for graphene and are strongly dependent on the state of 
the surface. The carrier concentration and correspondently the work 
function of graphene are exceptionally sensitive to substrate 
(intrinsic) and environmental (extrinsic) gating due to its two- 
dimensional nature. Being representative of the surface state, the 
work function can be strongly affected by such factors as traces of 
gas contamination, temperature and, in particular, humidity. For 
example, by conducting experiments in the controllable humidity 
(outside of the scope of the present paper) we found that standard 
(—10-15%) variations corresponding to typical day-to-day change 
of the lab humidity, led to a corresponding change of Ay CPD 12LG ~ 
25-30 mV. Moreover, it was reported that 1LG and 2LG have dif- 
ferent adsorption energies for gases, such as N0 2 22 and water 
vapour 15 , which leads to the differences in the doping levels for these 
two domains. Thus, without reproducing substrate and envir- 
onmental conditions, the SP and work functions of graphene 
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obtained in different experiments cannot be adequately compared. 
This places an even stronger emphasis on the necessity of standard- 
ization in the measurements of the SP and work function for each 
particular sample of interest and for given environmental conditions. 

Methods 

Sample preparation. Nominally monolayer epitaxial graphene was prepared by 
sublimation of Si and the subsequent graphene formation on the Si-terminated face of 
an on-axis 4H-SiC(0001) substrate at 2000°C and 1 bar argon gas pressure. Details of 
the growth and structural characterization are reported elsewhere 2 . The specific 
synthesis route has been developed to provide large areas of homogeneous single- 
layer graphene. The resulting material is n-doped, owing to charge transfer from the 
interfacial layer 23 ' 24 , with the measured electron concentration in the range n = 6-20 
X 10 11 cm" 2 and carrier mobility of ~ 3000 cm 2 V _1 s _1 at room temperature 25,26 . 

The epitaxial graphene device was fabricated by electron beam lithography 
(PMMA/MMA and ZEP520 resists), oxygen plasma etching and evaporation of Ti/ 
Au (5/100 nm) electrodes. Details of the sample fabrication are reported elsewhere 26 . 
The device comprises two crosses with a channel width of 4.8 urn, surrounded by 
1.6 um-wide trench etched down into the SiC substrate. The transport measurements 
were performed in air, at room temperature, in a dark environment. Details of the 
measurements are reported in Ref. 26. 

Standard lithography fabrication methods lead to a thin (1-2 nm) layer of a resist 
residue on top of the graphene. The residues can significantly affect the carrier density 
and even type, as shown in our previously published work in this area 27 . For instance, 
we show exposure of the resist residue to 250-nm wavelength UV light for 20 minutes 
leads p-doping of the graphene, where n h = 2.5 X 10 12 cm" 2 with 0 1LG = 4.68 eV. 
Subsequently, cleaning the residues using contact-mode AFM restored the n-type 
conduction of the graphene, where n e = 1.17 X 10 12 cm -2 with @hg = 4.35 eV. To 
avoid resist related doping, the device was cleaned by sweeping away the residual 
resist and partly atmospheric adsorbates from the surface using contact-mode AFM 
prior to imaging. In order to avoid permanent damage to the device, soft contact- 
mode cantilevers (Bruker) with a set point of —40 nN was used. 

SPM measurements. SP measurements were carried out in ambient environment at a 
controlled temperature of 18°C and humidity of —35%. The measurements were 
conducted on a Bruker Dimension Icon SPM. Doped silicon PFQNE-AL probes 
(Bruker) with a probe radius of —5 nm and a spring constant of —0.8 N/m were used 
for electrical measurements. Topography height images of the graphene device were 
recorded simultaneously with tapping phase and SP maps were compiled from either 
AM-KPFM, FM-KPFM or EFM phase shift. 

Electrostatic force microscopy. Electrostatic force microscopy (EFM) is performed 
as a dual-pass technique: first, the topography line profile is recorded in tapping 
mode, and then the line profile is traced at a set lift height above the surface. During 
the second lifted pass, the cantilever is mechanically oscillated at/ 0 , while a constant 
DC bias (V DC ) is applied, probing the probe-sample electrostatic forces, which 
depend on the probe-sample capacitance C and height z 28 : 

F DC =]- d ^V 2 , where V = V CPD + V DC + V induced . (1) 
2 az 

EFM is a purely DC technique. The electrostatic forces affect the amplitude, resonant 
frequency and phase of the cantilever oscillation. The EFM image is generated by 
recording the cantilever phase changes with a lock- in amplifier 5 

a QdF DC Q(d 2 C\ 2 

A(p= k^z- = TkW) v ^ (2) 

where k is the spring constant and Q is quality factor of the cantilever. The EFM 
technique operates on the force gradient (dF DC /dz) 29 , giving sharper contrast between 
areas of different electrical properties. Being confined to the probe apex, the force 
gradient decays much faster with distance than the force itself and, therefore it is less 
affected by the parasitic capacitance of the cantilever base. However, EFM provides 
only qualitative information on the electronic properties of sample surface, as the 
individual voltage components are not separated 29 . 

Electrostatic force spectroscopy. Electrostatic force spectroscopy (EFS) is performed 
at points of interest defined by EFM or other mapping techniques. Each measurement 
consists of oscillating the probe at/ 0 , while sweeping V probe and simultaneously 
recording Atp. The plots of Aq? as a function of V probe are parabolic, where the 
inflection point of the parabola is the point at which dF DC /dz is nullified, i.e. the force 
on the probe is zero (Figure lc). The inflection point is extracted post measurement 
and the resulting V probe at which dF DC /dz = 0 defines the surface potential. EFS 
spectroscopy was conducted along the center of the device channel, i.e. 200 
spectroscopy points were taken on the graphene channel along the marked line 
connecting the gold leads with the step of —300 nm between individual points. The 
step between individual points can be significantly decreased with the restricting 
factor being the lateral step resolution of the SPM system, which is generally limited 
by the diameter of the probe apex. Thus, the step wise change <20 nm in the surface 
potential at 1-2LG interface can be readily observed 30 . 



EFS can be used as a quantitative and accurate measure of the SP and Sample of a 
sample, if the probe is first calibrated against a sample of known 0. As EFS is not a 
scanning technique, probe degradation and the relevant work function change are 
negligible. EFS could be performed at every point of a two dimensional raster if time is 
not a constraint. 

Calibrated work function measurements of graphene were obtained with EFS by 
calibrating the work function of the probe against the known work function of gold 
electrodes, which was measured by ultraviolet photoemission spectroscopy (UPS), see 
Figure Id. 

Amplitude-modulated KPFM. The AM-KPFM, discussed here, is performed as a 
dual-pass technique; topography line profile is mapped with tapping mode AFM 
during the first pass, which is then traced at a set lift height above the surface 
performing the surface potential measurement (Figure la). During the second pass of 
AM-KPFM, the mechanical drive to the cantilever is disabled and an AC bias voltage 
(Vac = 2 V) is applied to the probe at the mechanical resonance f 0 of the cantilever. 
The V AC causes the cantilever to oscillate due to the attractive and repulsive 
electrostatic interaction (F es ) between the probe and the sample 31 

F es =- l - d ^- [(Vdc - V CPD ) + V AC sin (wtf (3) 
2 az 

where V DC is a DC bias voltage and V C pd is a contact potential difference between the 
probe and sample. A proportional-integral-derivative (PID) feedback loop monitors 
and minimizes the amplitude of the cantilever oscillations by applying a 
compensating V DC to the probe to cancel the probe- sample electrostatic forces, i.e. 
Vdc = V C pd is maintained at each pixel. The applied V DC is recorded at each point, 
providing a map of the SP. This conventional dual-pass KPFM is a well-established 
technique, widely used for quantitative probing of the surface potential of 
graphene 1 ' 3 ' 6 ' 32 ' 33 . Generally, the technique suffers from a poor lateral resolution, 
-50-70 nm, see e.g. Ref. 31. 

Frequency-modulated KPFM. FM-KPFM, discussed here, is a single-pass technique, 
which gives a greater degree of spatial resolution than AM-KPFM as it measures the 
force gradient (dF es /dz) 29 rather than the force acting on the entire cantilever. The 
probe- sample electrostatic forces affect the resonance frequency of the cantilever, 
where the amplitude (A) of the cantilever excitation at/ 0 — fmod depends on the 
electrostatic force in the system: 

A(f 0±/mod )=/ 0 (l-l§), (4) 

where k is the spring constant of the probe. The topography is determined with the 
tapping mode at the cantilever resonance, f 0 — 300 kHz. Simultaneously, a lower 
frequency (/mod ~ 2 kHz) AC voltage (V A c = 8 V) is applied to the cantilever. This 
modulation results in the appearance of side lobes in the cantilever oscillation 
spectrum at frequencies f 0 ± f mod (Figure lb). The FM-KPFM PID feedback loop 
minimizes the side lobes by applying a compensating V DC at each pixel. In a similar 
fashion to AM-KPFM, the probe- sample electrostatic forces are nullified when V DC = 
V C pd> therefore recording V DC and generating the SP map. However, in contrast to 
AM-KPFM, FM-KPFM typically requires stiffer, higher frequency cantilevers. FM- 
KPFM offers a higher spatial resolution of <20 nm as a result of force gradient 
localized to the probe apex and higher sensitivity to frequency shifts 29 . 

While AM-KPFM is usually performed as a dual-pass technique where first 
topography and then SP are measured along the same line in an alternating fashion, 
FM-KPFM is most often performed as a single-pass technique, where topography and 
potential are recorded simultaneously, thus improving the speed of image capture. 
However, it should be noted that being either single- or dual-pass is not a definition of 
the techniques, as other examples have been demonstrated previously 34 ' 35 . 
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